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ABSTRACT 

High-energy gamma-ray emission is theoretically expected to arise in tight binary star systems (■with high 
mass loss and high velocity ■winds), although the evidence of this relationship has proven to be elusive so 
far Here we present the first bounds on this putative emission from isolated Wolf-Rayet (WR) star binaries, 
WR 147 and WR 146, obtained from observations with the MAGIC telescope. 

Subject headings: gamma rays: observations 



1. INTRODUCTION 

WR stars represent an evolved stage of hot (Teff > 
20000 K), massive (Mzkus > 25 Af©) stars, and display 
some of the strongest sustained winds among galactic ob- 
jects: their terminal velocities may reach the range Voo > 
1000 — 5000 km/s. Perhaps with the exception of the short 
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lived luminous blue variables (LBVs), they have the high- 
est known mass loss rate M - 10^"*. ..10-^ Mo/yr of 
any stellar type. Thus, colliding winds of massive star bi- 
nary systems are considered as potential sites of non-thermal 
high-energy photon production, via leptonic and/or hadronic 
process after accele ration of prima ry particles in the col- 
Usion shock (e.g., Eichler & Usov I Il993l) . This possibil- 
ity is substantiated by the detection of (non-thermal) syn- 
chrotron radio-emission from the expected colliding wind lo- 
cation in some binaries (see below). Many models have 
been proposed to predict GeV to TeV emission from these 
binaries, w ith different levels of detail (e.g., among recent 



works see iBenaglia et al. II2OOII: iBenaelia & Romero 



iPittard & Dougherty I l2006t iReimer. Pohl & Reimer I 
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200i). 



Conceptual ly, the process would mimic the cases of L S 5039 
&iaronian et al. I BOOSbl). PSRB 1259-63 ([Aharonia n et al.l 
120061) or LS I H-6 1 303 ('Albert sLaLl EQOI, ^08d|ia), partic- 
ularly if they result in pulsar-driven 7-ray binaries in which 
7-rays may arise from a shock region produced by the inter- 
action of the winds of the two components. 

In a recent paper, the High Energy Stereoscopic array 
(H.E.S.S.) collaboration reported the discovery of very high 
energy (VHE) 7-ray emis sion coincident with the young stel- 
lar cluster Westerlund2 (lAharonian et al. Il2007h . The High 
Energy Gamma Ray Astronomy (HEGRA) and the Major At- 
mospheric Gamma Imaging Cherenkov (MAGIC) telescope 
detected the source TeV J2032-I-4130 and suggested a possible 
connection ■with the Cygnus OB2 cluster (Aharonian e t alH 
I2OO2I I2005at [Albert et al. I I2008ah . Theoretically, the rela- 

We recall that these binaries have orbital modulated TeV emission, and 
are formed by a massive star with a strong wind and a compact object, which 
only in the case of PSR B 1250-63 is known to be a pulsar. 
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tionship between stellar associations and high-energy emis- 
sion has been put forward by many authors, since early-type 
stellar associations have long been proposed as cosmic-ray 
acceleration sites and also as providers of target material for 
cosmic-ray interacti ons (see for instance the recen t papers by 
Parizot et al. 20 04: Torres. Domingo-Santamarfa & Rqmerql 
2004; Bednarek ll2005l:lDomingo-Santamarfa & Torres |2006|) 
and references therein. 

In the case of Cygnus OB2, the VHE emission is supposed 
to occur at a region displaced from the center of the asso- 
ciation, where detailed multiwavelength studies reve aled an 
overdensity of hot OB stars, although not WRs ( But t et al. I 
I2003L 120061) . In the case of Westerlund2, the stellar- cluster 
contains at least a dozen early-type O-stars, and two remark- 
able WR stars, WR20a and WR20b. In particular, WR20a 
was recently es tablished to be a binary (iRauw et al. 1 120041 : 
iBonanos et al. 1 12004). Based on the orbital period, the min- 
imum masses were found to be around (83 ± 5) Mcr^ an d 
(82 ± 5) Mq for the binai-y components (iRauw et al. Il2005h . 
what certainly qualifies it among the most massive binary sys- 
tems in our Galaxy. No significant flux variability or orbital 
periodicity was found in the corresponding data sets, neither 
for the Cygnus OB2 nor for the Westerlund 2 regions. Unless 
such an orbital period is detected in future datasets, it would 
be very difficult or impossible for the current generation of in- 
struments to distinguish whether the radiation observed from 
these associations is coming from a isolated binary system or 
rather is generated as a collective effect of the whole cluster 
Thus, a direct measurement of single WR binary systems, to 
explore whether they are able to produce high-energy 7-ray 
emission in an isolated condition, is worth pursuing. After 
briefly motivating the scenario for 7-ray emission from iso- 
lated binary systems, we present the results of MAGIC obser- 
vations of two such systems. 

2. HIGH-ENERGY 7-RAYS FROM WR BINARIES: 
CANDIDATES, OBSERVATIONS AND RESULTS 

For the present investigation, and apart from other technical 
considerations like a favorable declination, we looked for can- 
didates with non-thermal emission, indicating the presence of 
relativistic electrons, and for which the geometry of the col- 
liding wind region is established. The two selected systems 
-WR 146 and WR 147- have been resolved using the Very 
Large Array (VLA) and the Multi-Element Radio Linked In- 
terferometer (MERLIN) in two sources each and at least for 
WR 147, where the most detailed model is currently available, 
the system was predicted to be a powerful MAGIC source for 
most of the orbital period ( Reimer, Pohl & Reimer 2006). 

For the angular resolution of the MAGIC telescope (^ 
0.1°) the colliding wind zone will not appear to be spatially 
resolved, presenting individual colliding wind binary systems 
as point-source candidates at the 7-ray sky. Thus we have 
searched for point sources in the direction of these two bina- 
ries. 

2.1. WR147 

WR147 (see e.g.. ISetia Gunawan et aTl l2001l and refer- 
ences therein), among the closest and brightest systems that 
show non-thermal radio emission in the cm band, is composed 
of a WN8(h) plus a BO. 5 V star with a bolometric luminos- 
ity of Lboi = 5 X IO'^Lq and effective temperature Tcff = 
28500 K, i.e. thermal photon energy of about er ~ 6.6 eV. At 
a distance of 650 pc the implied binary separation is estimated 
to be 417 AU. The mass loss rates (Mwr = 2.5x IQ-^Mg/yr, 
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Fig. 1. — Inverse Compton (IC) spectra of WR 147 for orbital phases 0, 
0.25, 0.5 and 0.75 , neglecting any eccentricity of the system and assuming 
% = 90deg, fromlReimer, Pohl & Reimer] j200^ . 77 pair production ab- 
sorbs not more than <0.3% (> 50GeV) and <18% (> 100 GeV) of the 
produced flux at orbital phases 0.25 and 0.5, respectively. No absorption 
takes place at phase 0. MAGIC upper limits on this system are marked. 



Mob = 4 X lO^^AfQ/yr) and wind velocities (wwr — 
95Gkm/s, vqb = 8GGkm/s) place the stagnation point at 
6.6 X 10^"' cm. This is in fact in agreement with MERLIN 
observations, which show a northern non-thermal component 
and a southern thermal one with a separation of ( 575± 15) mas 
(IChur chwell eTaLlflOgl IWilUams et al. Ill997h . This radio 
morp hology and spect rum support a colliding wind scenario 
dWilliams et al. |[T997!) . whose collision region has also been 
detected by the Chandra X-ray telescope (Pittard et al. 2002). 

The non-thermal flux component can be well fitted by a 
power law with spectral index a = —0.43. Neither the 
eccentricity nor the inclination of the system are known 
due to the very long orbital period (1350 yr, as derived by 
[Setia Gunawan et al. 2001). Reimer, Pohl & Reimer (2006| 
have provided a detailed modeling of the high-energy 7-ray 
emission expected from WR 147. The expected fluxes are 
shown in Figure [T] together with MAGIC upper limits for 
which we give further details below. 

WR 147 was observed with the MAGIC telescope between 
1 1 August and 10 September 2007, for a total of 30.3 hours of 
good data (after quality cuts removing bad weather runs). The 
zenith angle of the observations ranged between 10° and 30°, 
being sensitive to gamma-rays in the energy range between 
about 80 GeV and lOTeV. The observations were carried out 
in the false-source track (wobble) mode ( IFomin et alJl 19941) . 
with two directions at 24' distance east-west of the source 
direction. The analysis of the d ata was performed w ith the 
MAGIC standard analysis chain dAlbert et al. Il2008ch . which 
combines Hillas image parameters J3y means of a Random 
Forest algorithm dAlbert et al. 1 2008b) for signal/background 
discrimination and energy estimation. The training of the al- 
gorithm is done by means of contemporary data from the Crab 
Nebula observations and Monte Carlo simulated gamma-ray 
events. Since February 2007, MAGIC signal digitization 
has been upgraded to 2GSample/s Flash Analog-to-Digital 
Coverters (FAD Cs), and timing param eters are used during 
the data analysis dTescaro et al. Il2007h . This results in an im- 
provement of the flux sensitivity from 2.5% to 1 .6% (at a flux 
peak energy of 280 GeV) of the Crab Nebula flux in 50 hours 
of observations. 

Searches of gamma-rays from WR 147 have been per- 
formed for three different energy cuts, namely: above 80 GeV, 
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TABLE 1 
WR 147 OBSERVATION RESULTS" 



Energy 


-^excess 


S 


U.L. 




[GeV] 


[evts] 


[t] 


[evts (cm^^ 


^)] 


>80 


-196±175 


-1.1 


150(1.1 X 10" 


-11) 


>200 


-92±89 


-1.0 


84 (3.1 X IQ- 


-12) 


>600 


-20±24 


-0.8 


28 (7.3 X IQ- 


-13) 



" From left to right: energy range, number of excess 
event s, statistical significance of the excess iLi&Ma] 
[T983) . and signal upper lim it for the different obser- 
vation nights. Upper limits jRolke. Lopez & Conrad I 
I2005h are 95% confidence level (CL) and are quoted 
in number of events and (between brackets) in 
photon flux units assuming a Crab-like spectrum 
iAlbert et al. I2008d) . 



TABLE 2 
WR 146 OBSERVATION RESULTS'' 



Energy 


^ excess 


S 


U.L. 




[GeV] 


[evts] 


[<t] 


[evts (cm~^ s~ 


^)] 


>80 


264±97 


2.7 


840 ( 3.5 X IQ- 


-11) 


>200 


133±67 


2.5 


487 ( 7.7 X 10" 


-12) 


>600 


-21 ±26 


-0.8 


46(5.6 X 10" 


-13) 



See explanation in Table[T] 

above 200 GeV and above 600 GeV. In all cases the number of 
signal candidate events found are compatible with statistical 
fluctuations of the expected background. The obtained upper 
limits are shown in Table [T] and Figure [T] and correspond to 
1.5%, 1.4% and 1.7% of the Crab Nebula flux for the three 
considered energy bins, respectively. 

2.2. WR146 

WR 146 is a similar system: a WC6H-08 coUiding- 
wind binary system also presenting thermal emission from 
the stellar winds of the two stars, and bright non- 
the rmal emission from the wind-collision region (e.g., 
see Dougherty et al. "1996'; Dougherty. Williams. & Pollaco 



|2000; O'Connor et al. 2005). The period is estimated to be 
^ 300 yr (Dou gherty et al. 1 119961) and the estimates of the 
distance to the system diff er from 0.75 kpc to 1.7kpc (see 
ISetia Gunawan et al. Il200ll and references therein for other 
system parameters). 

WR 146 is located ^ 0.7° away from the unidentified 
VHE 7-ray source TeV J2032H-4130 and was observed with 
MAGI C within the obser vation program devoted to this 
source jAlbert et al. Il2008ah . albeit with a reduced sensitivity. 
The total effective exposure, which accounts for the loss of 
sensitivity of off-axis observatio ns and camera ill umination 
during moonlight observations (lAlbert et al. I I2007D . is 44.5 



hours, obtained between 2005 and 2007. At MAGIC site, 
WR 146 culminates at 14° and the observations were carried 
out at zenith angles between 14° and 44° (see ' Albert et al. I 
l2008ai for details). The data analysis follows the standard 
MAGIC analysis chain. Since most of the data are acquired 
with 300 MHz FADCs, image timing parameters are not used 
in this analysis. 

The result of the searches of gamma-rays from WR 146 for 
three different energy cuts (above 80 GeV, above 200 GeV 
and >600 GeV) are shown in Table |2l As for the case of 
WR 147, all measured signal candidates are consistent with 
background fluctuations and the upper limits (corresponding 
to 5.0%, 3.5% and 1.2% of the Crab Nebula flux) are pre- 
sented. At the lower energy bins we see a positive number of 
excesses at the ~ 2.5cr level. However, the present data are 
too scarce to establish if this comes from a gamma-ray sig- 
nal or background fluctuations. Future observations will shed 
light on this issue. 

3. CONCLUDING REMARKS 

Our search for VHE 7-ray emission from two archetypi- 
cal cases of WR binaries produced the first bounds on such 
systems. These bounds constrain theoretical models (or, as- 
suming correctness of the models, their internal parameters, 
such as the -unknown- orbital phases of the systems). The 
establishment of WR binaries as VHE 7-ray sources is yet 
pending. 

The case for WR 147 as a potential 7-ray source 
for MAGIC was theoretically established before, as 
shown in the correspo nding curves of Figure [1] from 
iReimer. Pohl & Reimer I (12006 ). for most of the orbital 
phases. The validity of this model is baselined on an assumed 
ensemble of orbital parameters, which are still unknown for 
this system. For instance, ignorance of its current phase as 
well as of its eccentricity and inclination makes a direct rul- 
ing out of this model impossible, although that the presented 
scenario could nominally survive only for phases close to 0, 
defined where the line of sight encounters first with the WN8 
and then the B0.5V star. The MAGIC observations show that 
irrespective of phase, GLAST should see a flux cutoff well 
within its range of detectability in the tens of GeV regime, if 
it is able to detect the stars at all. 
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curves depicted in Figure 1. We also would like to thank the 
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in La Palma. The support of the German BMBF and MPG, the 
Italian INFN and Spanish CICYT is gratefully acknowledged. 
This work was also supported by ETH Research Grant TH 
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